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Abstract 

The Multi Module Model for Ultra-Relativistic Heavy 
Ion Collisions at RHIC and LHC energies is presented. It 
uses the Effective String Rope Model for the calculation of 
the initial stages of the reaction; the output of this model 

*Talk given at the New Trend in High-Energy Physics, Yalta, Crimea, 
Ukraine, September 22-29, 2001. 



1 



is used as the initial state for the subsequent one-fluid cal- 
culations. It is shown that such an initial state leads to 
the creation of the third flow component. The hydrody- 
namical evolution of the energy density distribution is also 
presented. 

A realistic and detailed description of an energetic heavy ion re- 
action requires a Multi Module Model, wherein the different stages 
of the reaction are each described with a suitable theoretical ap- 
proach. It is important that these Modules are coupled to each 
other correctly: on the interface, which is a three dimensional 
hyper-surface in space-time, all conservation laws should be satis- 
fied and entropy should not decrease. These matching conditions 
were worked out and studied for the correct matching at FO in 
detail in Refs. 0. 

Fluid dynamical models are widely used to describe ultra- 
relativistic heavy ion collisions. Their advantage is that one can 
vary flexibly the Equation of State (EoS) of the matter and test its 
consequences on the reaction dynamics and the outcome. For ex- 
ample, the only models which may handle the supercooled QGP 
are hydrodynamical models with a corresponding EoS. In ener- 
getic collisions of large heavy ions, especially if a Quark-Gluon 
Plasma (QGP) is formed in the collision, one-fluid dynamics is a 
valid and good description for the intermediate stages of the re- 
action. Here, interactions are strong and frequent, so that other 
models, (e.g. transport models, string models, etc., that assume 
binary collisions, with free propagation of constituents between 
collisions) have limited validity. 

On the other hand, the initial and final, Freeze-Out (FO), 
stages of the reaction are outside the domain of applicability of 
the fluid dynamical model. After hadronization and FO, the mat- 
ter is already dilute and can be described well with kinetic models. 
The initial stages are the most problematic. Non of the theoreti- 
cal models on the physics market can unambiguously describe the 
initial stages (see for example 0, |^ for discussion). 

Thus, since we can not unambiguously describe the initial 
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stages of relativistic heavy ion collisions in microscopic models, we 
have to use some phenomenological models. There are two phe- 
nomenological models, which are most frequently discussed in the 
literature - the Landau and the Bjorken scenarios. These models 
describe two extreme cases - complete stopping or complete trans- 
parency, respectively. Nowadays, experiments have entered the 
region where the Bjorken model is expected to be applicable, but 
there is no clear and unambiguous confirmation that experiments 
follow this scenario. The preliminary experimental results from 
RHIC do not show transparency - most particle multiplicities (ex- 
cept maybe the most central collisions) do not show a dip in the 
rapidity spectra, but rather a plateau around mid rapidity |^ , 
which is a sign of strong stopping. Furthermore, very strong ellip- 
tic flow has been measured, which shows a clear peak around mid 
rapidity |^ . To build such a strong elliptic flow, strong stopping 
and momentum equilibration are required. Also the p/p ratio at 
mid rapidity measured at RHIC ||TD|, ^ (preliminary) is still far 
from 1, which tells us that the middle region is not baryonfree. 

Our goal is to build a Multi Module Model for ultra-relativistic 
heavy ion collisions valid for RHIC and LHC energies, and maybe 
for the most energetic SPS collisions. In the present work only 
the first step is done - a Two Module Model ||, |^ - but a very 
important one, since the Effective String Rope Model (ESRM) 
[0, 1^ has been developed for the most problematic module - the 
module describing the initial stages of collisions. 

We describe the initial state with the ESRM; then the hydro- 
dynamical evolution with a QGP EoS is used for the intermediate 
state. The hydrodymanical calculations are performed with the 
Los Alamos Particle-in-Cell (PIC) one fluid code O, The hy- 



dro evolution stops at the FO hypersurface. We present a version 
of the code assuming that FO (simultaneous chemical and ther- 
mal) happens on the simplified toothed hypersurface, where it's 
normal vector, da^, is parallel to the flow velocity for every cell. 
On average this hypersurface approximates the constant time hy- 
persurface. Therefore the flow velocity does not change during the 
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FO process, and the calculations can be done in local rest frame 
of the matter. Such a surface is also completely timelike; this al- 
lows us to avoid the problems discussed in Refs. |l|]. The more 
advanced description of the FO process is planned to be separated 
into the Third Module. 

The EoS presently used in the code was used and presented in 



Ref. 1 13]. This phenomenological phase transition model 



A) takes the phenomenological EoS for hadronic matter in a simple 
form, which nevertheless allows one to check different parameter- 
izations discussed in the literature; 

B) uses the Bag model EoS to describe QGP; 

C) creates a complete EoS, containing pure phases and a region 
where they coexist, by the Maxwell construction. 

The reader may notice that the model is still in a formative 
stage and a lot of further work remains necessary. We are not yet 
ready to present the quantitative calculations to be compared with 
data. Nevertheless we present some preliminary results. First of 
all, the flow component, fi, as a function of rapidity, y, is pre- 
sented in Fig. to show that our expectation that the model will 
generate a third flow component (Refs. 0) is reasonable [0, ^. 
Our initial state generated by the ESRM indeed produces a strong 
antifiow in semi-central collisions. The directed vi component ap- 
pears to be very small, as expected for RHIC energies (for example 
[^]). The peaks in Vi (actually v^) around y = ±0.3 look very high, 
but this is typical for calculations without thermal smearing [|1^ 



Including thermal smearing will lead to smaller and wider peaks. 

We also present in Fig. ^ the evolution of the energy density 
distribution, e{r,t), during the hydrodynamic expansion (notice 
that t = corresponds to the moment then we start the hydro 
description, i.e. 3.5 fm/c after the collision). We can see the 
development of the small but very dense, shock-like domains at 
the intermediate stages. This effect, at least to some extent, is due 
to numerical artifacts. The PIC method has a well-known defect 



called 'ringing' or marker-clustering [16|. It is called this way 



because early calculations in two dimensions resulted in density 
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distributions that looked similar to the vibration patterns in sound 
vibrations - areas of higher and lower densities Advanced 



algorithms that mitigate this effect are being implemented in the 
code. 

Another reason which may cause the development of the nu- 
merical artifacts is a step-like in beam direction initial energy den- 
sity distribution (output of ESRM): it has a jump from e{x, y) = 
const inside the matter to in the outside vacuum. In order to 
avoid (or at least to suppress) the effect of this we are planning to 
smooth over initial energy density distribution, for example by a 



Gaussian shape, as proposed in Ref . [0 . Thus it is plausible that 
the development of the shock-like domains in the simulation is not 
a physical effect, although we can not rule out other possibilities. 

To evaluate the observables a third, FO module is going to be 
attached to the model. Here the widely used Cooper-Frye model 
should be essentially improved and modified. First of all conser- 
vation laws and the requirement of increasing entropy should be 
enforced in the module. In addition, particularly for FO across 
space-like hyper surf aces, realistic, non-equilibrium post FO phase 
space distributions will have to be used to avoid negative contri- 
butions occuring in the naive use of the Cooper-Frye model [|l], |^ . 
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Figure 1: Preliminary results of Two Module Model calculations 
for vi vs rapidity, y. Au-|-Au collision at Eq = 65 GeV/nud for 
different impact parameters. Initial state was generated with A = 
0.065, T = 3.5 fm/c (see g | for details). 
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Figure 2: Preliminary results of Two Module Model calculations 
for energy density distribution in the reaction plane. Au+Au col- 
lision at = 65 GeV/nucl, b = 0.1 ■ 2i?^„ for different moments 
of time; initial state was generated with A = 0.08, r = 3.5 fm/c 
(see H, i for details). g 



